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Abstract 

The interaction of the low lying pseudo-scalar mesons with the ground state baryons in the charm sector is 
studied within a coupled channel approach using a t-channel vector-exchange driving force. The amplitudes 
describing the scattering of the pseudo-scalar mesons off the ground-state baryons are obtained by solving 
the Lippmann-Schwinger equation. We analyze in detail the effects of going beyond the t — approximation. 
Our model predicts the dynamical generation of several open charmed baryon resonances in different isospin 
and strangeness channels, some of which can be clearly identified with recently observed states. 
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I. INTRODUCTION 



A very active topic of research in hadron physics concerns the study and characterization of 
resonances, in order to establish whether they qualify as genuine qq or qqq states or, alternatively, 
they behave more as hadron molecules generated dynamically. A series of pioneer works 

naaa 

[5,(3], based on a t-channel vector meson exchange force, already predicted a wealth of s- wave baryon 
resonances generated by coupled channel dynamics with effective hadronic degrees of freedom rather 



than quarks and gluons. The earlier approaches have been readapted in the 
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ast decade to the 
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22, 



23, 
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modern language of chiral lagrangians [3, |8|, [9 
and many resonances in the light SU(3) sector, which cannot be described properly by quark 



models 25J unless substantial meson-baryon components are included [26j |. have been identified 
with dynamical states generated from the interactions of mesons of the pseudoscalar 0~ octet with 
the l/2 + ground state baryons. Some consequences of these studies, such as the two-pole nature 



of the A(1405), are confirmed through the analyses 



29, 



27, 



30, 



281 ] of different experimental reactions 



311 ] . Note that the basic structure of a molecular- type baryon is quite different from that 



implied by the quark models, even when they include the dressing with meson-baryon components, 
as in the 3 Pq formalism, where the qqq Hilbert s pac e is coupled to the meson-baryon Hilbert space 



through the creation of a qq pair 



32, 



33 



34 



35j, [36], [37]] • In the former case the degrees of freedom 



are purely hadronic and the resonances must be seen as pseudo-bound states of two hadrons, while 
in the later case the essential component of a baryon is still of three quark nature. A goal in 
hadron physics research is to distinguish between both pictures from thorough analyses of as many 
properties of the hadron as possible, such as the mass, width, magnetic moment, form-factors, 
etc. . . 

In recent years it was demonstrated that besides the s-wave baryon resonances many more 
states can be generated dynamically. Baryon resonances with J p = 3/2- were studied based 
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39, 
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411 ] . D-wave 



on the leading order chiral lagrangian with the decuplet 3/2 + fields 

jaryon resonances were also generated dynamically with vector meson degrees of freedom in Refs. 



42 



43, 



44 



45 , [47(- Another promising line of research is the recent interpretation of low 
l/2 + resonances as molecular states of two pseudoscalar mesons and one baryon 
(5^]. All these results support the so-called hadrogenesis conjecture, formulated 
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ymgj 



48, 



49 



50 



51 



a few years ago by Lutz and Kolomeitsev, according to which resonances not belongin g to the 



large- iV c ground state of QCD are generated by coupled-channel dynamics [la . |42j, 



53, 



54, 



55]. 



The study of charmed hadrons is receiving an increased attention thanks to the efforts of a 
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series of collaborations, both at lepton colliders (CLEO, BELLE, BaBar) and hadron facilities 
(CDF@Fermilab, PHENIX and STAR@RHIC, FAIR@GSI). The new results confirm with better 



statistics previously seen charmed states but are alsogivmg rise to t 



of new hadrons 
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60, 



61 
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63 



64, 
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66 



67 



68, 



re discovery of a large amount 



69l | . Coupled-channel unitary 



schemes have been recently extended to include the charm degree of freedom and have been applied 
to the description of open and hidden charm mesons, with the o 
a straightforward interpretation as meson molecules 
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72, 



aservation that some states admit 
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Similar methods have 



been used for describing baryons with charm, motivated in part by a clear parallelism between the 



behavior of the A (1405 



one 
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80 



in the C = 0, S = -1 sector with the A c (2595) in the C = 1 and S = 
?e consistent with the spin-flavour Heavy Quark Symmetry that 



81|- To 



develops in this heavy sector 



82 



83 



841 ] . the vector mesons and J = 3/2 + baryons have recently 



jeen included in the basis of meson-baryon states, emp 



oying a static spin-flavour SU(8) scheme 



similar to that developed in the light sector [43j, |44|. Treating the D and D* mesons on an 
equal footing basis has lead to the observation that some of the dynamical generated states, such as 
the A c (2595), have mostly a D*N composition rather than the DN molecular nature. In any case, 
the observation that some of the dynamically generated charmed hadrons can be readily identified 



with observed resonances, such as the J p = 1/2 A 



baryons, sustains the hadrogenesis conjecture 



m, 



42 



2595) or the J p = 3/2~ A c (2625) charmed 



53, 



54 



Apart from the different basis of states included in the models, a common feature of all the 
previous works is the use of an interaction based on the i-channel exchange of vector mesons, as the 
driving force for the s— wave scattering of pseudo-scalar mesons off ground state baryons. The limit 
t — * is applied, leading to a vector type WT zero-range interaction. This procedure is justified 
for on-shell meson-baryon transitions, MB — * M'B 1 , which are diagonal (M'B' = MB) and hence 
the value of t is small as long as one is not too far from threshold. It also holds for non-diagonal 
amplitudes {M'B' ^ MB) that show a moderate difference of masses between the initial and final 
mesons and baryons involved, as is the case of meson-baryon scattering within the light SU(3) 
world. However, the coupled channel dynamics in the heavy sector finds also charm-exchange 
processes for which the difference of masses between the external mesons are comparable with the 
mass of the charmed vector meson being exchanged. This clearly signals the breakdown of the 
zero-range approximation which is no longer reliable for these non-diagonal transitions. While one 
may still argue that many of the dynamically generated states are triggered by a single dominant 
meson-baryon interaction component and, hence, their energy can be well estimated by the pole 
position of an uncoupled calculation involving only diagonal amplitudes, the corresponding width 
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will however be determined by non-diagonal amplitudes and will therefore depend on whether the 
t = approximation is implemented or not. Moreover, it is well known that some resonances 
owe their origin to a particularly strong coupling between different channels, hence involving non- 
diagonal transitions, in which case the t = approximation is not at all appropriate for these 
states. 

In the present work we study the charmed baryon resonances obtained dynamically from 
the interaction of the low lying pseudo-scalar mesons with the ground state baryons within a 
coupled channel approach, using the full t-dependence of the t-channel vector-exchange driving 
term, instead of the t = approximation. We incorporate the t-dependence within a general 
four-dimensional integration scheme, which we reduce to a three-dimensional equation of the 
"Lippmann-Schwinger" type, and analyze in detail the effects of going beyond the t = ap- 
proximation within this scheme. 

The paper is organized as follows. In Sec. [II] we present the formalism, showing the details of 
the kernel employed and the equation used to obtain the scattering amplitudes. Our results for 
the properties of the baryon resonances with charm in various strangeness and isospin cases are 
shown in Sec. IIIIl where we also compare with the t = results. A summary of our conclusions is 
presented in Sec. IIVI 



II. FORMALISM 



78] , we identify a t-channel exchange of vector 



Following the original work of Hofmann and Lutz 
mesons as the driving force for the s-wave scattering. In their original model, Hofmann and Lutz 
exploited the universal vector meson coupling hypothesis. They considered the t-channel exchange 
of vector mesons between pseudoscalar mesons in 16-plet and baryons in 20-plet representations in 
such a way to respect chiral symmetry for the light meson sector. The Weinberg-Tomozawa type 
interaction is recovered in the zero range limit, (i.e., t — > and see Eqs. (6) and (8) of Ref. [3] for 
details). The scattering kernel has the form 

ve[i6] v v / v 

(1) 

where the sum runs over all vector mesons of the SU(4) 16-plet, (p, K* , K* , u>, 4>, D* , D*, 
D*, D*, J/ty), my is the mass of the exchanged vector meson, g is the universal vector meson 
coupling constant, ki, qi, kj and qj are the four momenta of the incoming and outgoing baryon and 



4 



meson, and the coefficients C^'.y denote the strength of the interaction in the different sectors 
(I, S, C = 1) and channels The value of g = 6.6 reproduces the decay width of the rho meson 



Assuming that \t/my\ « 1, the t-dependence of the scattering kernel can be neglected, giving 

rise to the zero-range approximation. The s-wave projection of the scattering kernel under such 

approximation is easily obtained, and in the center-of-mass frame it takes the form 

,(I,S,C) 
'ij;V 

4 , m? 



2 (J^- 

V^\h,q u k v q 3 ) = -N 9 - Y, U\k i \)+E(\k i \)+u(\k j \)+E(\k j \)-M i -M j 



ve[ie] 



"v 



mh — mf \ 
-- L Z^^- M i) \ > ( 2 ) 

where rrii,mj, Mj, Mj are the masses of the incoming and outgoing mesons and baryons, and 
tji(\ki\),ujj(\kj\), Ei(\ki\), Ej(\kj\) their corresponding energies, which have been taken to be their 
on-shell values. The factor N = [(E(\ki\) + Mi){E(\kj\) + M j )/(4M i M i )] 1 / 2 comes from the nor- 
malization of the Dirac spinors. Note that for physical (fully on-shell) transition amplitudes one 
has Ui(\ki\) + Ei{\ki\) = ujj{\kj\) + Ej(\kj\) = y/s, in which case one recovers the familiar expression 
of the commonly used Weinberg-Tomozawa interaction. The last term in Eq. ([2]) is usually ignored 
in most works, to be more consistent with the t = approximation applied in the denominator 
of the meson-exchange propagator, as noted in Ref. 



only minor corrections 



811 ] . In any case, its consideration introduces 



In order to illustrate the validity of this approach we show in Fig. Q] the value of t/rriy for 
cos# = — 1 as a function of y/s, where my is the mass of a representative meson exchanged, which 
we take to be the p meson mass for diagonal transitions and the D* meson mass for charm exchange 
ones. The range of energies goes roughly between the 7rE c and DN thresholds, thereby covering the 
region of the J p = l/2~ resonance A c (2595), which is a prime example of a dynamically generated 



78 



79 



80, 



811 ] . It also expands beyond 



open charmed baryon state in various approaches [76|,l77 
the DN threshold for about 300 MeV in order to explore the energy region that will be relevant 
in future studies of the D-meson self-energy in a nuclear medium. As one can see, the value of 
t/rriy is only close to zero for diagonal transitions around their corresponding energy threshold, 
and its size can be comparable to one at the energies of interest. For the non-diagonal 7rS c — > DN 
transition, t/rriy never goes to zero and acquires values of the order of 0.5. 
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FIG. 1: Dependence on the center-of-mass energy *Js of the four momentum tranfer t/niy for cosO = —1 
and for different transition amplitudes. 

The results of Fig. Q] point clearly to the need for exploring the effects of going beyond the t = 
approximation, an attempt that is taken in the present work by considering the full i-dependence 
of the scattering kernel given by Eq. (pQ). Performing the s-wave projection, we obtain the following 
analytic expression: 



2 r 2/3 ah-Pa. (a + b 



ve[i6] 



H 77T — In 



b b 2 \a-b 



(3) 



with a, b, a and (5 being 



mf + m| - 2^(1^1)^(1^1) - 



& — 2 1 1 1 kj | 

2 2 
jYl . — TTli- 
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Uly 
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2 9 \ 

777, . — 777, - \ 

*-(n j -(|A; i |)-ni(|fc i |)+M J --Af i ) , (4) 



m 



where we have defined f2(|fc|) = cj(|A;|) + E(\k\). 

The i = expression of the s-wave interaction is recovered by expanding the logarithm of Eq. ([3]) 
in the limit b/a — ► up to the linear term in b/a, and setting a = —m v . As one can infer from the 



values of t/rriy displayed in Fig. [fl keeping the t-dependence in the denominator of the kernel [see 
Eq. (H|)] will decrease the strength of the diagonal transitions, since in this case t < 0. Therefore, 
in order to reproduce a given resonance found in local models, the present approach will in general 
need to compensate the lack of strength with a higher value of the cut-off momentum used to 
regularize the loop integrals. On the other hand, non-diagonal amplitudes, responsible mainly for 
the decay width of the dynamically generated states, will be enhanced since they are characterized 
by a positive time-like t value due to the large mass difference between the mesons and baryons 
involved in the transition. As a consequence, our resonances will be wider than those found in the 
local models. 

Retaining the f-dependence in the kernel implies additional analytical structures which 
prevent us from obtaining the scattering amplitudes by solving the Bethe-Salpeter (BS) equation 
using on-shell amplitudes only. Instead, we incorporate the t-dependence within a more gen- 
eral four-dimensional integration scheme, which we reduce to a three-dimensional equation of the 
"Lippmann-Schwinger" type. To this end, for a given value of the total scattering energy, we will 
evaluate the transition potential between any arbitrary pair of relative momenta within the cut-off 
value, keeping the dependence on the momentum transfer in the corresponding matrix element. 
This procedure follows the same spirit of the usual meson-exchange models of the NN interaction 
also applied to meson-baryon scattering models 



89 



90|. We note that keeping the full t de- 
pendence in the exchanged meson propagator also implies that retardation effects are implemented 



as done in the three-dimensional reduction of the Bethe-Salpeter equation of Ref. 



911 ]. which is 



different than the prescription based on time-dependent perturbation theory [88l |. Both choices 
differ in the way off-shell effects are implemented and their differences in the NN sector show up 
especially in the contributions of the two-meson exchange box diagrams [88f] not included in the 
present model. In any case, observables can be matched to experimental data with either choice 
of retardation effects by selecting appropriate values of the renormalization parameters. 

We note that Eqs. jl]) to (|3|) assume infinitely lived (zero width) exchanged vector mesons, while 
some of them have large widths due to their strong decay into a pair of mesons, like the p meson 
contributing to diagonal channels. We have checked, however, that the value of t is never larger 
than the square of the minimum energy required for the meson to decay, namely (2m T ) 2 in the 
case of p exchange or (m n + mx) 2 for K* . In other words, the mesons being exchanged in this 
problem are largely off-shell and they will be treated as stable particles. 

Once the scattering kernel has been constructed, we can obtain the T-matrices which describe 
the scattering of the pseudo-scalar meson fields off the baryon fields by solving the well known 
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Lippmann-Schwinger equation 

(5) 

where 

J&*& (VS,k) = =r+ , • (6) 

2E m (\k\)io m ( y \k\) yfi-E m (\k\) -io m (\k\) +irj 

We have introduced a dipole-type form factor 

F(\k\) = ( A \ ] (7) 
\A 2 + \k\ 2 J 

to regularize the integral. This form is typically adopted in studies of hadron-hadron interactions 
within the scheme of Lippmann-Schwinger-type equations in the light flavour sector 



The 



value of the cut-off A is a free parameter of our model. Given the limited amount of data for 
charmed baryon resonances, and in order to simplify the analysis, the cut-off A is adjusted to the 
position of a well known J p = l/2~ state in a particular isospin and strangeness sector, and the 
same value is used for the other sectors explored in this work. We will also investigate the effect 
of a gaussian-type form factor, as well as the dependence of our results on the value of the cut-off 
employed. 

Note that the approaches based on the Bethe-Salpeter equation solved with on-shell amplitudes 
ignore the off-shellness (momentum dependence) of the kernel and scattering amplitude in the loop 
function. Actually, the on-shell factorization can only be justified within the WT form of the 
potential which is obtained after applying the t — * limit. 

Another aspect worth commenting is the different extrapolation of the kernel at subthreshold 
energies. While the kernel of the on-shell BS approaches depends on -y/s in the form given by 
Eq. j2|) upon replacing tUi(\ki\) + Ei(\ki\) + Uj(\kj\) + Ej(\kj\) by the potential used in the 
tri-dimensional Lippmann-Schwinger equation depends essentially on the incoming and outgoing 
tri-momenta. These are always taken real quantities in our approach, hence the sum of the four 
single particle energies is always larger than the sum of the two meson-baryon thresholds involved 
in the transition. An explicit dependence on yfs is only implemented in the meson-baryon interme- 
diate propagator of the Lippmann-Schwinger equation that determines the scattering amplitudes. 
Therefore, at subthreshold energies, the factor in the numerator of the kernel used in on-shell BS 
approaches is smaller than in the present work. This in part compensates the enhancement in 
diagonal transitions associated to the t — > limit. In any case, the free parameters of the model 



(cut-off values) can finally be conveniently fine-tuned to adjust the energy position of a well known 
resonance, as is usually done for the case of the A c (2595) appearing about 200 MeV below the 
threshold of the channel DN to which it couples very strongly. In this respect, the differences 
between the present work and the on-shell BS approaches will be made more evident in properties 
tied to non-diagonal transitions which are reduced in the t — > limit, as is the case of the resonance 
widths. 

In order to associate a given enhancement of the scattering amplitude to a resonance, we look 
for a characteristic pole in the unphysical sheet of the complex energy plane. Our prescription of 
unphysical sheet is such that, whenever the real part of the complex energy crosses a meson-baryon 
threshold cut, the sign of the on-shell momentum is changed for this channel and for the already 
opened ones, as described in detail in Ref. [4[. Once a pole zr is found, its value determines 
the Breit-Wigner mass (M = Hezn) and width (r = 2Imz/j) of the resonance, as seen from 
real energies, if the pole is not too far from the real axis. The couplings of the resonances into 
the meson-baryon components of a given sector are obtained from the residues of the scattering 
amplitude since, close to the pole, it can be parameterized in the form: 

T ;ff (k i ,k j>Z ) = ^-. ( 8) 

Note that, as it stands, the value of the coupling constants of Eq. ([8]) depend on the particular 
momentum values chosen in the evaluation of the T-matrix element. Since we are only interested 
on the size of the couplings relative to the various channels, we take the prescription of evaluating 
them for the case ki = kj = 0. 

III. RESULTS AND DISCUSSION 

One of the main interests of this paper is to study the effects of going beyond the t = 
approximation in the kernel, as has been customarily done. We have already anticipated in the 
previous section that, in the case of coupled channel problems involving light and heavy flavored 
particles, this approximation is not always justified. 

All the possible sectors with charm C = 1 that can be built from the s-wave scattering of 
pseudoscalar mesons with J p = l/2 + baryons are shown in Table HI together with the corresponding 
meson-baryon coupled channels. In this work, we will first study the cases in which some resonance 
with either J p = l/2~ or unknown spin-parity has already been observed. This includes the sectors 
with isospin, strangeness quantum numbers (I,S) = (0,0), (1,0) and (1/2,-1), corresponding 
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(i-3,i) 




(0,-2,1) 


KE C KE' C DE r)tt c r/fi c D s tt cc r/ c Sl c 


(1,-2,1) 


Trn c KE C KE' C DE C 


(1,-1,1) 


ttE c ttE' c KA c KZ c DA V E c DZ V E' c KVL c D s E V 'E c tj'E' c Vc E c D s E cc Dtt cc i h E' c 


(1,-1,1) 


ttE c ttE' c KY, c £>E 


(0,0,1) 


ttE c DN v A c KE c KE' c D s A 77' A c ?7 C A C DE cc 


(1,0,1) 


ttA c ttE c DN KE C r)Z c KE' c D s Z rj"S c DE CC 7? C E C 


(2,0,1) 


7rE c 


(5,1,1) 


KAc D S N KY, C 


(1,1,1) 





TABLE I: Coupled-channel meson-baryon states with charm C = 1 and all possible combinations of isospin, 
strangeness (I,S). 

respectively to A c , E c and H c states, the experimental information of which is gathered in Table ILTl 
We will next explore the sector (J, S) = (0, —2) of the Q c states which so far has no experimental 
evidences for J p = l/2~ states. Finally, we will comment on the (/, S) = (1/2, 1) sector that can 
only be realized with the presence of 5 quarks. 



Resonance [MeV] 


I(J P ) 


Width [MeV] 


Decay modes 


Status 


A c (2595) + 


o(|") 


3.6+J 


A C 7T7T, E C 7T 


*** 


A c (2765)+ or E c (2765) 


?(? ? ) 


- 50 


A C 7T7T 




A c (2940) + 


0(? ? ) 


17±S 


iVL>, E c tt 


*** 


E c (2800) 


1(? ? ) 


75lf 7 (E++), 62+«°(S+), 61±?«(E°) 


A C 7T 


*** 


S c (2790) 


KD 


< 15(5+), < 12(5°) 






5 C (2930) 


?(? ? ) 


36 ± 13 


A c if 




S c (2980) 


m 


26 ±7(5+), 20 ±7(5°) 


A c i?7T, E C K 




5 C (3055) 


?(? ? ) 


17 ± 13 


A c i?7T, E e K 


** 


5 C (3080) 




5.8 ± 1.0(5+), 5.6 ±2.2(5°) 


A c Ktt,Z c K,Z* c K 


*** 


S c (3123) 


?(? ? ) 


4±4 


E*X 





TABLE II: Masses, widths, decay modes and status of experimental charmed baryon resonances with J p = 
l/2~ or unknown. 
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A. A c resonances (I=0,S=0,C=1) sector 



In this sector there exists a three-star narrow J = 1 /2 resonance, the A c (2595), which has 



been extensively studied in various works 



76 



78, 



3, 3 



8l[ . We start by comparing in Fig. [2] the 



results obtained using our non-local kernel with those taking the limit t — > 0. We represent the 
imaginary part of the scattering amplitude of the elastic process DN — > DN, as a function of sfs 
for zero incoming and outgoing relative momentum values. We can see that, by adjusting the cut- 
off value conveniently, both models of the kernel can generate this state dynamically. However, the 
zero range approximation needs a cut-off value of A = 553 MeV while the finite range interaction 
requires a substantially larger value of A = 903 MeV. 

_ 1.0 



> 

CD 



^ 0.5 

z 

Q 

A 

z 

Q 



0.0 





1 

— t*0 


■ 


-- t=0 _ 







2590 



2595 
s 1/2 [MeV] 



2600 



FIG. 2: Imaginary part of the scattering amplitude of the elastic process DN — > DN in the (I, S, C) = 
(0, 0, 1) sector as a function of y/s for the finite range interaction (solid line) and the zero range approximation 
(dashed line). The incoming and outgoing relative momenta ki and kj have been taken equal to 0. 

This is easily understood from the fact that the DN — ► DN diagonal matrix elements of the 
non-local potential, largely responsible for generating the resonance, are smaller in magnitude 
than those of the local one. The large difference between the cut-off momentum values is just a 
reflection of the importance of the non-local terms in this problem. Once the A c (2595) resonance 
is conveniently located to its experimental position by both prescriptions, there remain substantial 
differences in its width. The local potential produces a very narrow resonance, of width 0.15 MeV, 
while the resonance generated by the finite range potential has a width of 0.5 MeV, closer to the 
empirical value of 3.6 + 2.0 — 1.3 MeV. Again, this is due to the different magnitude of the non- 
diagonal matrix elements DN — > 7rS c , which are larger in the finite range approach. We note that 
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our model does not consider the three-body decay channel A c 7T7t which already represents almost 
one third of the decay events 



86l | . We also observe that the results obtained with the low cut-off 
in the approximation t = for the A c (2595) agree (mass, width and couplings) with former studies 
of meson-baryon resonances of Hofmann and Lutz in Ref. [78] and Garcia-Recio et al, in Ref. 851 ] . 

Our search of resonances in this sector produces two states that are listed in Table [TTTl together 
with their widths and couplings to the various meson baryon states. We immediately see that the 
A c (2595) is basically a DN state which couples very weakly to its only possible decaying channel 
7rS c , thereby explaining its narrowness. We obtain an even narrower resonance at 2805 MeV, which 
is a K S c bound system, a state also found around the same energy in Ref. |78j and Ref. 



. Note 



that this resonance couples non-negligibly to DN and, was its location be moved upward in energy 
by 20-30 MeV with a slight change of the cut-off parameter, it could explain part of the structures 
seen below 2.85 GeV in the D°p invariant mass spectrum measured by the BaBar collaboration 
68|. Table [TTTl also shows the results obtained with the local t = model. We remark that, in spite 
of the fact that the second resonance appears at a higher energy, 2827 MeV, its width is narrower 
than in the finite-range model, confirming the trend observed for the A c (2595). 



(I = 0,S = 0,C= 1) 



A [MeV] 


903 (t ^ 0) 


553 (t = 0) 


M [MeV] 


2595 


2805 


2595 


2827 


T [MeV] 


0.5 


0.01 


0.15 


0.006 


Couplings |<7j 




tt£ c (2591) 


0.44 


0.001 


0.21 


0.002 


LW(2806) 


16.34 


0.23 


18 


0.03 


?7A C (2832) 


0.58 


1.92 


0.22 


1.60 


K5 C (2963) 


0.71 


3.76 


0.23 


3.22 


KS' C (3Q7Q) 


0.32 


0.004 


0.01 


0.13 


D S A(3085) 


8.16 


0.18 


8.54 


0.07 


r/A c (3243) 


0.96 


0.01 


0.54 


0.003 


?7 C A C (5265) 


2.83 


0.02 


1.74 


0.007 


5S CC (5307) 


0.07 


0.96 


0.03 


0.48 



TABLE III: Masses, widths and couplings of the resonances in the (I, S, C) = (0, 0, 1) sector, for the non-local 
(t 7^ 0) and local (t = 0) models. 



We represent in Fig. [3] the modulus square of the coherent sum of all transition amplitudes 
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FIG. 3: Modulus square of the coherent sum of all transition amplitudes going to any of the possible final 
meson-baryon decaying channels, as a function of y/s, for the (J, S, C) — (0, 0, 1) sector. Our results are 
split into two panels, covering different energy regions and having different energy scales, to better visualize 
the properties of each state. The incoming and outgoing relative momenta fej and kj have been taken equal 
to 0. 



going to a final meson-baryon state to which the resonances can decay, namely: 

I Cm, b ,T M i B i^ mb (\/s) | 2 , 



(9) 



M'B' 



where the values of the coefficients C^, B , would depend on the specific reaction used to excite 
the resonance and, in this graphical example, have all been taken to one. The amplitudes have 
been calculated for zero incoming and outgoing relative momentum values. To compare to actual 
experiments, one should be using the appropriate excitation coefficients as well as transition matrix 
elements going to the on-shell final momentum corresponding to the value of y/s. Therefore, the 
results in Fig. [3] and similar ones throughout the paper should be considered mereley illustrative. 
The representation is split into the various energy regions where the resonances appear. In general, 
a resonance couples dominantly to a given channel and the value of the maximum of Eq. © is 
basically proportional to the modulus squared of the product of the resonance couplings to the 
dominant and decaying channels, gu'B' an d §mb respectively, and inversely proportional to the 
resonance width T. Note that, instead of adjusting the vertical axis of Fig. [3}o to the maximum 
associated to the narrow resonance at 2805 MeV, we have scaled it down to better visualize the 
enhancement at 3069 MeV, right below the KB! C threshold. This enhancement becomes a resonance 
if we increase the cut-off value slightl y, w ith properties that are similar to the state found around 



the same energy by the local models 



78, 



85|. 
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B. E c resonances: (I=1,S=0,C=1) sector 



{1 = 1, 5 = 0,(7 = 1) 


A [MeV] 


903 (t ^ 0) 


553 (f = 0) 


M [MeV] 


2551 


2804 


2585 


2804 


T [MeV] 


0.15 


5 


0.005 


0.63 


Couplings \gi\ 


ttA c (2424) 


0.06 


0.27 


0.002 


0.04 


ttE c (2591) 


4.00 


0.16 


2.15 


0.04 


£W(2806) 


1.25 


2.10 


0.38 


1.70 


i^S c (2963) 


0.04 


0.20 


0.003 


0.06 


?7S C (2999) 


0.79 


0.11 


0.44 


0.03 


#5^(3070) 


2.30 


0.14 


1.55 


0.04 


£> s £(3162) 


0.62 


1.79 


0.17 


1.37 


7/E c (3410) 


0.04 


0.19 


0.006 


0.09 


5~ cc (5307) 


0.91 


0.15 


0.30 


0.02 


t? c £ c (5432) 


0.13 


0.55 


0.02 


0.27 



TABLE IV: Masses, widths and couplings of the resonances in the (I, S, (7) = (1, 0, 1) sector, for the non-local 
(t 0) and local (t = 0) models. 

Using the same cut-off that reproduces the A c (2595) in the (J, S, C) = (0, 0, 1) sector, we predict 
two narrow resonances at 2551 and 2804 MeV that appear right below the thresholds of the channels 
to which they couple mostly, namely ttT, c and DN, respectively, as can be seen from Table IIVI 
We also show the results obtained with the t = model and the same cut-off value of A = 553 
MeV adjusted to the position of the A c (2595) in the isoscalar channel. As expected, the resonance 
at 2804 MeV, which couples mostly to DN states, barely changes its position from that of the 
non-local model. Indeed, since the A c (2595) also couples mostly to DN, the cut-off adjustment of 
the local model using this resonance as reference has essentially left the DN amplitude intact in 
this energy region, thereby generating similar "DN-type" bound states as in the non-local model 
in the various isospin sectors. Note, however, that the width of the resonance at 2804 MeV is an 
order of magnitude smaller in the local model. The lower energy resonance, coupling mostly to 7rS c 
states, appears at somewhat larger energies in the local model and, in spite of the gain in phase 
space, its width is substantially reduced. This comparison confirms the trend already observed in 
the case of the 1 = sector. The resonances of the local model appear at similar energies as those 
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of the non-local approach but their widths are much smaller. Since this will be a general behavior 
in all sectors, we will only show, in the remaining sections, the results obtained with the non-local 
approach developed in this work. 
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FIG. 4: The same as Fig. E] but for the (I, S, C) = (1,0, 1) sector. 



The sum of amplitudes squared is represented in Fig. HI Since the E c (2551) resonance couples 
negligibly to its only allowed strong decaying channel, ttA c , it is seen in Fig. 0^ as a narrow peak 
which makes it easy to miss given the energy resolution of the meson-baryon invariant masses 
built up in the present experiments. The couplings to the different meson-baryon states of the 
S c resonance at 2804 MeV, visualized in Fig. 0b, allow one to identify it with the state found in 
Ref. [3] at a substantially lower energy, 2680 MeV, using a subtraction method to regularize the 



loops, as well as that found in Ref. [8l|] around 2750 MeV, using a cut-off method which preserves 
isospin symmetry in the regularization scheme. Our result is obviously closer to that of the latter 
work. 

n 

The Belle Collaboration reported recently [62] an isotriplet of excited charmed baryons, 
S c (2800), decaying into A+7r~, A+7r° and A+7r + pairs and having a width of around 60 MeV 
with more than 50% error. Although this resonance has been tentatively assumed to decay to A c 7r 
pairs in d-wave and its spin parity estimated to be J p = 3/2~ , actual angular distributions have 
not been measured and the fits to A c tt spectra can not rule out s-wave type decays. Hence, our state 
at 2805 MeV could be easily identified with the S c (2800) resonance, provided three-body decay 
mechanisms not accounted for in our model, could explain the large width observed experimentally. 
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C. S c resonances: (1=1/2, S=-1,C=1) sector 



The results in this sector are presented in Table IVl and Fig. We obtain two pure bound states 
at 2515 and 2549 MeV, respectively, which are placed less than one pion mass above the H c member 
of the J p = l/2 + ground state antitriplet, (A+, S+, H^), and below the mass of the 'H! c member 
of the J p = l/2 + sextet, (E°, £+, H' c °, H' c + , ft®). This implies that these bound states would 
decay electromagnetically through the emission of H c 7 pairs and may have been detected at photon 
energies of about 50 and 80 MeV in the experiment where the r! c was observed 63J. Although 



no apparent signals are reported, we note that the limited statistical significance of the spectra 
measured in [63}] prevents one from ruling out the existence of these bound states. Moreover, their 
production rate would also be somewhat inhibited by the fact that they are predominantly 5 quark- 
component states. Note also that, lowering gradually the value of the cut-off to somewhat below 
700 MeV, the first state at 2515 MeV eventually becomes resonance, but quite narrow due to its 
weak coupling to the first channel, whereas the second state at 2549 MeV, which is a 7rH c molecule, 
rapidly becomes so wide that it would be difficult to be distinguished from the background. 

In addition, our model gives three more resonances above the 7rH c threshold and below 3 GeV, 
placed at 2733, 2840 and 2977 MeV. The local model of Ref. [3], based on on-shell amplitudes, 
also obtains three resonances in this energy region, located in general at somewhat lower masses 
and showing a different order of appearance, as can be inferred from the values of their couplings 
to the different meson-baryon components. More specifically, the lowest resonance in the local 
model appearing at 2691 MeV and coupling strongly to DT, should be identified with our middle 
resonance at 2840 MeV. The next two resonances appear quite close in the scheme of Ref. [3], at 
2793 MeV and 2806 MeV, coupling mostly to KT, C and DA, respectively, while in our case they 
are further apart from each other, at 2733 and 2977 MeV. The crossing in the ordering of states 
is another consequence of the different values of the transition potential amplitudes used in both 
coupled-channel schemes. 

Recently, several S c states have been observed by the CLEO 



6J, Belle |65( and BaBar 



66] 



Collaborations, out of which the three star possible candidates to be identified with one of our 
states are at 2790, 2980 and 3080 MeV [see Table HI]. A change of the cut-off value within a 
reasonable range could bring any of our two lower mass resonances to agree in position with the 
H c (2790) but the width would turn out to be twice wider than the observed one in the case of the 
lower mass state. The H c (2980) could be easily associated to either one of the two higher mass 
states found here. However, the experimental analysis of Ref. jfjfj concludes that the H c (2980) + 
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(1 = 1/2,5= -1,C = 1) 


M [MeV] 


2515 


2549 


2733 


2840 


2977 


r [MeV] 


0. 


0. 


34 


0.58 


4 


Couplings |fifj| 


ttS c (2609) 


0.65 


4.47 


0.05 


0.06 


0.31 


7^(2715) 


4.84 


0.76 


1.77 


0.01 


0.22 


#A C (2779) 


0.48 


3.21 


0.19 


0.10 


0.19 


#£ c (2946) 


6.90 


1.01 


7.37 


0.93 


0.16 


L>A(2985) 


1.03 


0.30 


0.96 


1.54 


2.95 


r?S c (3018) 


0.13 


1.04 


0.13 


0.18 


0.10 


L>£(3062) 


2.91 


0.89 


3.64 


8.82 


1.74 


^(3124) 


4.04 


0.59 


3.47 


0.46 


0.07 


*TQ C (3192) 


4.40 


0.68 


1.47 


0.07 


0.21 


D S H(3288) 


1.74 


0.15 


0.92 


4.71 


2.51 


t?'H c (3428) 


0.16 


0.05 


0.15 


0.53 


0.06 


^(3534) 


0.01 


0.03 


0.09 


0.11 


0.28 


5 S S CC (5408) 


0.02 


0.01 


0.25 


0.35 


0.71 


5fi cc (5429) 


1.10 


0.53 


1.18 


0.24 


0.03 


?7 C S C (5450) 


1.17 


1.10 


0.55 


0.04 


0.03 


77^(5556) 


0.04 


0.08 


0.26 


0.33 


0.82 



TABLE V: Masses, widths and couplings of the resonances in the (J, S, C) = (5,-1, 1) sector 
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FIG. 5: The same as Fig. El but for the (I, S, C) = (1/2,-1, 1) sector. 
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state decays in about 50% of the cases into T,^ + K pairs, which makes our state at 2840 MeV, 
showing a stronger coupling to KH C , the most likely candidate to be associated to the S c (2980). 



D. fi c resonances: (I=0,S=-2,C=1) sector 



In this sector we predict the existence of a bound state at 2959 MeV, near the lowest threshold, 
and two resonances placed at 2966 and 3117 MeV as can be seen in Table IVII and Fig. [6J The 
possible bound state could be detected through the decay into f2 c 7 states with photons of = 260 
MeV in the center-of-mass frame. The resonance placed at 2966 MeV and seen in Fig. [Oh, is very 
narrow (r = 1.1 MeV) according to the low coupling of the resonance to the only channel in which 
it can decay (Kz> c ) and the little available phase space. The resonance at 3117 MeV with a width 
of r = 16 MeV, seen in Fig. [6b, is a DH molecule that can decay into K3 C and K3' c states. 



(7 = 0,6' = -2, (7=1) 


M [MeV] 


2959 


2966 


3117 


T [MeV] 


0. 


1.1 


16 


Couplings \gi\ 


KS C (2964) 


1.36 


0.43 


0.51 


ifS' c (3070) 


2.04 


4.49 


0.27 


DS(3189) 


2.03 


1.68 


5.34 


77fi c (3246) 


1.67 


3.69 


0.24 


T7'S1 C (3656) 


0.10 


0.07 


0.35 


D s n cc (5528) 


0.17 


1.17 


0.19 


?/ c ft c (5678) 


0.28 


0.21 


1.03 



TABLE VI: Masses, widths and couplings of the resonances in the (J, S, C) = (0, —2, 1) sector 



The work of Ref. 78J also finds three states but placed at lower energies, 2839, 2928 and 2953 
MeV, which follows the trend observed for other sectors. The pattern of couplings to the various 
meson-baryon states also differs a little owing to the different interaction model used. The highest 



energy resonance in Ref. [78j, coup 
one, while the lowest one in Ref. 
resonance at higher energy. 



ing strongly to Krl c and r]{l c , would correspond to our middle 



781 ]. coupling strongly to DH, would be the equivalent to our 
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FIG. 6: The same as Fig. [3] but for the (I, S, C) = (0, -2, 1) sector. 
E. Resonances of five quarks 

We have also analyzed the sectors corresponding to resonances that cannot be realized with 
only three quarks and, therefore, their existence would be signaling the presence of pure five quark 
states. Note that the possible pentaquark-type systems predicted by the present model would be 
color singlet states built up from combinations of color singlet qq with color singlet qqq components. 
States with a different composition, such as a combination of color octet qq with color octet qqq 
clusters, can not be generated by our meson-baryon scattering approach. 



(7=1/2,5=1,(7 = 1) 


A [MeV] 


903 


1200 


1400 


M [MeV] 


2946 


2941 


2924 


r [MeV] 


0.93 


5 


12 


Couplings \gi\ 


iiTA c (2779) 


0.002 


0.04 


0.10 


D S N(2908) 


0.03 


0.84 


1.68 


if£ c (2946) 


0.07 


1.79 


3.59 



TABLE VII: Masses, widths and couplings of the resonances in the (I, S, C) — (1/2, 1, 1) sector for different 
cut-off values: 903, 1200 and 1400 MeV. 

Out of the three possible sectors, namely (I,S,C) = (2,0,1), (|,1,1) and (|,1, 1), we only 
find hints of a possible resonance in the case I = 1/2, S = 1,C = 1, where we see a cuspdike 
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FIG. 7: The same as Fig. [3] but for the (I,S,C) = (1/2,1,1) sector and for three different cut-off values: 
903, 1200 and 1400 MeV. 

structure placed at the threshold of the KTi c channel to which the state couples more strongly. 
This behavior is shown by the solid line in Fig. [7| and by the first column of coupling constants 
displayed in Table IVIIl which have been obtained using our nominal cut-off value of 903 MeV. 



According to the mechanism discussed in Ref. 



441 ]. the coupling constants should vanish if the 



resonance was placed right at the KT, C threshold, which explains the smallness of their values. In 
order to see whether the cusp structure would eventually become a clear resonance with a slight 
change of parameters, we also display in Fig. [7] and in Table PVTll our results with two other values 
of the cut-off, 1200 and 1400 MeV. One can clearly see that the cusp structure becomes a more 
bound and wider resonance as the cut-off value increases, while the coupling constants become 
larger. 

F. Dependence on model parameters 

We finalize this work by exploring the dependence of our results on the shape and size of the 
form factor employed, which are ingredients of the model that are not constrained by symmetry 
arguments. 

First, replacing the dipole-type form factor by a gaussian form: 

_ P 

F(\ k |) = e ™$ , 

we are able to adjust the position of the A c (2595) with a gaussian cut-off value of A 5 = 543 MeV. 
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The corresponding width is exactly the same as that found for the dipole-type form factor. When 
exploring the other sectors, the gaussian form factor gives rise to the same resonances, some of 
them slightly displaced by at most 50 MeV from the position found with the dipole-type form 
factor, but having essentially the same width. 

Retaining the dipolar form factor, we next explore the effects of varying the value of the cut-off 
A, a parameter that is not constrained by symmetry arguments, between 600 and 1200 MeV, that 
is, up to 300 MeV below and above the nominal value of 903 MeV used in this work. This variation 
produces changes in the positions and widths of the resonances within certain ranges, the general 
trend of which are summarized in the following points: 

1. A resonance that lies far below -by 50 to 200 MeV- the meson-baryon threshold to which it 
couples more strongly may change its position by an amount comparable with the variation 
of the cut-off value. The larger the cut-off the more bound the resonance becomes. 

2. The width of the resonance only changes appreciably for cut-off values that move the reso- 
nance above the threshold of a meson-baryon channel to which the resonance couples signif- 
icantly. 

3. Weakly bound resonances change their positions more moderately, at most by 10 MeV for 
changes of cut-off values within 100 MeV. In this case, the width tends to decrease as the 
resonance becomes less bound because of the distortions induced by moving closer to the 
threshold, a phenomenon also known as Flatte effect 921 ] . 



Having explored the systematics to the cut-off changes, we finally summarize in Table IVTlTI the 
states which, taking an appropriate cut-off value within the range explored, could be identified 
with a well established resonance of J p = l/2~ or unknown spin-parity. 



IV. CONCLUSIONS 



We have studied charmed baryon resonances obtained from a coupled channels unitary approach 
using a t-channel vector-exchange driving force. 

To the best of our knowledge, all previous models of dynamically generated baryon resonances 
in the charm sector rely on a local zero-range interaction, which is obtained by neglecting the four- 
momentum transfer t in front of the mass of the exchanged vector meson squared, m v . However, 
we have illustrated that the value t/m v is not at all negligible in the heavy sector, especially for 
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(I,S) 


A [MeV] 


Theory 

Mass [MeV] main channel Width [MeV] 


Experiment 
Mass [MeV] Width [MeV] Status 


(0,0) 
A c 


903 


2595 £W(2806) 0.5 


2595.4 ±0.6 3.6^3° *** 


(1,0) 

E c 


1100 


2792 ZW(2806) 16 


2801^ 751?? (£++) 
2792+^ &2+% (El) *** 
28021 4 6ll^ (E°) 




814 
980 

655 
960 


2790 ifE c (2946) 55 
2790 L>E(3062) 0.5 

2970 L>E(3062) 1.2 
2970 L>A(2980) 5.1 


2789.1 ±3.2 < 15 (Si) *** 
2791.8 ±3.3 < 12 (£°) 

2971.4 ± 3.3 26 ±7 (S+) *** 
2968.0 ± 2.6 20 ±7 (S°) 



TABLE VIII: Masses, widths and main coupled channel of states that can be identified with well established 
resonances in various sectors. 

charm-exchange amplitudes which produce a large value of the four-momentum transfer due to the 
large difference bewteen the masses of the mesons involved in the transition. 

We have analyzed in detail the effects of going beyond the t = approximation and, taking the 
1 = 0, C = l,J = l/2 sector of the well established J = l/2~ A c (2595) resonance as reference, we 
find that the experimental data is better reproduced by the non-local model which also requires a 
more reasonable cut-off regularization value of 903 MeV. 

Compared to the local models based on on-shell amplitudes, our approach obtains basically the 
same amount of resonances in all sectors but appearing, in general, at somewhat larger energies 
because the diagonal amplitudes, largely responsible for generating the bound states, are smaller in 
magnitude. An essential finding of this work is that our non-local approach produces much wider 
resonances because of the larger value of the non-diagonal amplitudes when t ^ 0. 

Varying the cut-off parameter within a reasonable range, we are able to locate some of our 
states at the energy position of a measured resonance in the same sector. In particular, we suggest 
the identification of the A c (2595), S c (2800), H c (2790) and ~ c (2980) as dynamically generated 
resonances having J p = l/2~ . In general, the widths of the states produced by our model are 
smaller than the experimentally observed ones, since we do not account for three-particle decay 
channels. 

We find a possible resonance in the sector with quantum numbers (I, S, C) = 1, 1) that can 
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only be realized by the consideration of a minimum of five quarks. The cusp-like structure observed 
at the threshold of the KT, C channel for a cut-off value of 903 MeV, becomes a more bound and 
wider clear resonance as the cut-off value increases. 

This is the first exploratory study of the effects tied to the non-locality of the meson-baryon 
interaction in the charm sector. We have considered meson-baryon coupled states built up from 
the J p = 0~ mesons and the ground state J p = l/2 + baryons. However, Heavy Quark Symmetry 
demands that the states containing heavy vector mesons are treated on equal footing due to the 
similarity of the masses with the heavy pseudoscalar ones. This will be addressed in a future work. 
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